Condition for the Existence of Complex Modes in a Trapped Bose— Einstein 
Condensate with a Highly Quantized Vortex 
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We consider a trapped Bose-Einstein condensate (BEC) with a highly quantized vortex. For the 
BEC with a doubly, triply or quadruply quantized vortex, the numerical calculations have shown 
that the Bogoliubov-de Gennes equations, which describe the fluctuation of the condensate, have 
complex eigenvalues. In this paper, we obtain the analytic expression of the condition for the 
existence of complex modes, using the method developed by Rossignoli and Kowalski [R. Rossignoli 
and A. M. Kowalski, Phys. Rev. A 72, 032101 (2005)] for the small coupling constant. To derive 
it, we make the two-mode approximation. With the derived analytic formula, we can identify the 
quantum number of the complex modes for each winding number of the vortex. Our result is 
consistent with those obtained by the numerical calculation in the case that the winding number 
is two, three or four. We prove that the complex modes always exist when the condensate has a 
highly quantized vortex. 

PACS numbers: 03.75.Kk,05.30. Jp,03.75.Lm,03.70.+k, 



I. INTRODUCTION 

The Bose-Einstein condensates (BECs) of neutral 
atomic gases were realized in 1995 0, 0, 0, and several 
kinds of vortices inside the BECs have been observed: 
the singly quantized vortex Q , the vortex lattice 0, @ , 
and the doubly quantized vortex which was created by 
the topological phase engineering Q. In particular, an 
interesting phenomenon was observed: a doubly quan- 
tized vortex decayed into two single quantized vortices 
spontaneously @- 

The theoretical investigations on the instability of the 
highly quantized vortex of the neutral atomic BEC were 
made by several authors 0, 0> E3 ■ They solved the 
Bogoliubov-de Gennes (BdG) equations [Tj, [Tj, [l5| nu- 
merically, and found that the equations have complex 
eigenvalues when a condensate has a vortex with the 
winding number two 0, j three 01 or f° ur G2 ■ The 
complex eigenvalues cause the blowup or d amp ing of the 
c-number condensate fluctuations lid 

[HGj. This in- 
stability of the condensate, caused by the complex eigen- 
values of the BdG equations, is called "dynamical insta- 
bility" . But it is still unknown whether a highly quan- 
tized vortex with an arbitrary high winding number al- 
ways brings complex eigenvalues and which eigenvalues 
become complex. It should also stressed that the rela- 
tion between the "dynamical instability" in theoretical 
concept and the observed decay of the doubly quantized 
vortex is not elucidated fully and still under investigation 
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As for the "complex modes", there is known another 
treatment, developed by Rossignoli and Kowalski [Isj . 
We refer to it as the RK method. In this method the 
quantum Hamiltonian of the quadratic form of creation 
and annihilation operators is considered, and the com- 
plex modes appear as a result of diagonalizing it with un- 
usual operators which are neither bosonic nor fermionic 
ones. Obviously the meaning of the complex modes in 
the RK method is different from one mentioned above 
(lfl ITU IT^ ] : they are quantum fluctuations in the for- 
mer, while they are c-number ones in the latter. The 
relation between the complex modes and the complex 
eigenvalues of the BdG equations has not been estab- 
lished. Further, we point out that the interpretation of 
the complex eigenvalues of the BdG equations in the case 
that they are treated as quantum fluctuations, as well as 
of the complex modes of the RK method, is not simple. 
We recently proposed a new approach of quantum field 
theory, using the eigenfunctions of the BdG equations in- 
cluding those belonging to the complex eigenvalues and 
considering the free quantum Hamiltonian consistently 
[lflj . although the content of this paper is not related to 
the approach directly. 

In this paper, we derive the conditions for the ex- 
istence of the complex modes in a trapped BEC with 
a highly quantized vortex analytically using the RK 
method within the two-mode approximation and small 
coupling expansion. Using the condition, we can answer 
the question above: does any highly quantized vortex 
bring the complex modes? The answer is yes, and we 
can identify partially which mode is the complex one. 
The condition is also useful for checking the results of 
the numerical calculation. 

This paper is organized as follows. In Sec. [H] the 
model Hamiltonian describing the weak interacting neu- 
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tral atoms in a harmonic trap is introduced. The c- 
number field whose square represents the distribution of 
the condensate is subject to the Gross-Pitaevskii (GP), 
and the quantum field is expanded in terms of a complete 
set. In Sec. lIIII we give a brief review of the RK method. 
In Sec. HVI first, we explain the two-mode approximation. 
Second, the small coupling expansion is introduced. Both 
of them play essential roles in our analysis. We derive 
the condition for the existence of the complex modes an- 
alytically within the two-mode approximation and small 
coupling expansion. Finally, we prove that the complex 
modes always exist when the condensate has a highly 
quantized vortex. Section is devoted to the summary. 



II. MODEL HAMILTONIAN 

We start with the following Hamiltonian to describe 
the neutral atoms trapped by a harmonic potential of a 
cylindrical symmetry, 



H 



d A x 



ft(x)(K + V{r, z)-(i)fax) 



+ 9 -fa\x)fa{x) 



where r = \/ x 2 + y 2 , x = (x,t) and 



K = V 2 , 

2M ' 

V(r, z) = ~M(a,ir 2 + lo 2 z z 2 ) 



(1) 

(2) 
(3) 



with the mass of a neutral atom M, the chemical poten- 
tial n and the coupling constant g. The bosonic field op- 
erator ip(x) obeys the canonical commutation relations, 

[fax), fa{x')]\ t=t , = S(x - x>), (4) 
[fax), fax')\\ t = t , = [fa(x), fa(x')]\ t ^ =0. (5) 

Let us divide the field operator ip(x) into a classical field 
£(x) and a quantum field <p(x) as 



fax) = £(x) + <p{x) 



(6) 



Here, the classical field £(x) is the order parameter repre- 
senting the condensate with a (highly) quantized vortex, 
which is assumed to be time-independent. Note that the 
function !;(x) is essentially complex due to the existence 
of the quantized vortex. At the tree level, £(x) satisfies 
the following GP equation: 

[K + V(r,z) - M + g |£(x) | 2 ] = . (7) 

The condensate particle number iV c is given by 



d 3 x\ax)\ 2 = n c 



(8) 



Substituting Eq. JBJ into Eq. QJ, one can rewrite the 
total Hamiltonian H as 

H = H + ffint , (9) 

where the free and interaction Hamiltonians are chosen, 
respectively, as 



Hn = 



dh 



ip f {x) {K + V(r, z) - fi} tp(x) 



+ y -{A\i{x)\ 2 0\xmx) 

+ e 2 (x)0 2 (x)+e(x)0 n (x)} 



(10) 



9 / d 3 x 



t{x)^ 2 {xmx)+e(xW(x)<p 2 (x) 



1 



(11) 



Assuming that a vortex is created along the z-axis and 
taking account of an axial symmetry along the vortex 
line, we introduce a real function /(r, z), defined as 



t{x) = J^e™ e f{r,z) 
V Ztt 



(12) 



where the integer k > 1 is a winding number of the vor- 
tex. The function f(r, z) is normalized as 

rdrdzf 2 (r,z) = 1 . (13) 

Substituting Eq. (|12"|) into Eq. J7J, one obtains 

, + V(r,z)-/i+^/ 2 (r,z) f(r,z)=0, (14) 

where 



2M \3r 2 + rdr r 2 + dz 2 ) ' ( j 
Using the bosonic creation and annihilation operators 
%i 1 m anc l ®n,t,m> we expand the bosonic field <£>{x), 

oo oo 

& x ) = J2 an,t,me i{ - l+K)e u nAm {r,z), (16) 

n,m— t— — oo 

where n (= 0, 1, 2, • • • ) is the principle quantum number 
(representing the radial nodes), t (= 0, ±1, ±2, • • • ) is the 
magnetic quantum number, m (= 0, 1, 2, • • • ) is the quan- 
tum number along the z axis. The function u n ^^ m {r, z) is 
a solution of the eigenequation with an eigenvalue e n ^ iTn , 



[K K+t + V(r, z)] u n ^ m (r, z) = (r,z), (17) 

and satisfies the following orthogonal and completeness 
conditions: 



2tt J rdrdzu n ^ m (r,z)u n >^ m '{r,z) = 5 n . n i5 m . m , , (18) 

V" u n ,e.m(r, z)u n .t. m (r', z) = - — 5(r - r')S(z - z') . 

(19) 



n,m—0 
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The eigenf unctions in Eq. 117(1 are expressed with the La- 
guerre polynomials L^(x) and the Hermite polynomials 

u nAm (r, z) = C nAm e-^ r \a^rf+ K \L^{alr 2 ) 
x e-i a > 2 H m (a z z) , (20) 



where 

C n ,l,r. 



tt 3 /2 V 2 m to! (\1 + k\ +n)l 



(21) 



a± = y / Mu± and a z — \[MuT z . The eigenvalues are 
given as 

e n ,e,m = w_l (2n + \l + k\ + 1) + u) z (m + ^ . (22) 

The free Hamiltonian is a sum of the sectors W l \ labeled 
by the absolute value of the magnetic quantum number 
j = |*|: 



Ha = & = H° + H 3 ' 

3=0 j' = l 



(23) 



Hereafter we use the symbols j as j = 0, 1, 2, • • • , and j' 
as j' = 1', 2', • • • , to make it explicit whether the special 
sector with j — is included or not in the arguments. 
The sectors H° and H 3 (j' > 1) are explicitly given 

as 



E 



n.n' ,m,m'—0 



«n m & n',0,m'A(n, °> m 5 n> i °> 



+ g {Sn,o,man',o,m'B(n, 0, to; n', 0, to') + h.c.} 

"(24) 

OO 

H 3 ' = ^2 ^n,j-,m & n',j\m'A{n,j\m\n',j',m') 

n,n f ,m,m f =0 

+ K_ rtm a n ',-j',m'A(n, -j', to; n', -j', to') 

+ ^{anj',man',-j\ m 'B(nJ',m;n',-j',m') + h.c.} 



1 



2 

+ h.c.} 



{an,-j', m a n 'j',m'B(n, -j',m; n',j',m') 



(25) 



with 



A(n,£,m;n ',£' ' ,m!) = {(e n ,i,m - fJ-)$n,n'Sm,m> 

+ 4K(n, £, to; n', £', to')}^ e , 

(26) 

B(n, £, to; < to') = 2K(n, £, to; n', m')S e ,-e , (27) 
where -ftT(n, £, to; n', to') is defined as 
K(n, £, to; n , £ , to') 



rdr dz / 2 (r, 2:)u„,f, m (r, zjv/.m'^ z ) • (28) 



III. ROSSIGNOLI-KOWALSKI METHOD 

In this section, we briefly review the RK method. In 
this method the "complex modes" appear in the diag- 
onalized Hamiltonian with unusual operators which are 
neither bosonic nor fermionic ones |18| . 

We begin with the general Hamiltonian for the BEC of 
the quadratic form of creation and annihilation operators 
a N and cln as 



N.N' ^ ^ 1 



+ ^{BNN'CL^a^, + h.c.) 



where 



= ^nz, 



n = 



(29) 

(30) 
(31) 



Here, the matrices A and B are hermite and symmetric, 
respectively. The symbols d and d^ represent the arrays 



A B 
B* A 1 



of components d^ and d 



t . 

AT" 



0>N=1 



'®N=1 
l N=2 



(32) 



(33) 



where N stands for the combination of the quantum num- 
bers n, I and to, and we treat N as a single integer 
for simplicity. The bosonic commutation relations of the 
original creation and annihilation operators djy and a N 
are written as 



where 



M = 



1 

-1 



(34) 



(35) 



Now, we consider a general linear canonical transfor- 
mation of the d and d^ 



z = wz' , 



Z' 



(36) 
(37) 
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where 



*JV=1 
1 N=2 



N=2 



(38) 



(39) 



We note that a' N is not necessarily the hermitian conju- 
gate of a' N , although a' N and a' N are required to satisfy 
the bosonic commutation relations as 



z'z' - (Z n Z n Y = M 



where 



Z' = 
T = 



z n r . 



1 

1 o 



As Z't = Z'W, the matrix W should satisfy 
WMW = M , 

where 

W = TW l T . 



(40) 

(41) 
(42) 

(43) 
(44) 



It is important to notice that in the case where only the 
real modes appear, the relation W = holds, and then 



N ■ 



The quadratic Hamiltonian i? qua d is rewritten as 

1 



H qu3 ,d — —ZHZ , 

h' = whw . 



(45) 
(46) 



Finding a representation in which H' is diagonal corre- 
sponds to solving an eigenvalue equation with the "met- 
ric" M, i.e., 



HW = MWMH' . 



(47) 



We regard this equation as an eigenvalue equation for a 
non-hcrmitian matrix TL: 



HW = WW , 
H = MH = 



A 

-B* 



B 

-A 1 



(48) 
(49) 



with the diagonal matrix W whose diagonal components 
are eigenvalues of H, denoted by xn- Then the quadratic 
Hamiltonian i? qua dratic can be expressed as 



H, 



quad 



N 



— I *> * 



(50) 



IV. CONDITION FOR THE EXISTENCE OF 
THE COMPLEX MODES 

Applying the RK method [lj§| reviewed in the pre- 
vious section to our Hamiltonian in Eqs. 1231) (I28|l . we 
now solve the eigenvalue equation (|48|l under the cir- 
cumstance to obtain the mode-decoupled form of our 
Hamiltonian. But the equation is too difficult to solve 
analytically in general. 

In this section, we develop the following two approx- 
imate schemes in Eq. ()48|l . First, we make a two-mode 
approximation on the Hamiltonian l|23|l . Second, we in- 
troduce the small coupling expansion. Both of them are 
helpful for the analytical study in spite of the nonlinearity 
of the basic equation for the condensate . And finally, 
we analytically derive the condition for the existence of 
the complex modes. 



A. Two-Mode Approximation 

In this subsection, we explain the two-mode approxi- 
mation which is used to obtain the analytic expression of 
the condition for the existence of the complex modes. 

First, we consider the j = part of the free Hamilto- 
nian 1)24(1. Within the two-mode approximation, H° is 
written in the matrix representation as 



ZqTLqZq 



(51) 



where 



Zn = 



7i 



7~io 

A 
Bo 



0,2 



V4/ 

B AoJ ' 

ei -n + 4K(l,l) 4K"(1,2) 

4#(2,1) £2-^ + 4^(2,2) 

(2K(1,1) 2K(1,2) 
\ 2K(2, 1) 2K{2,2) 



(52) 

(53) 
(54) 

(55) 

(56) 



Here, the mode 1 represents one labeled by the quan- 
tum numbers (n, £, m) — (n\, 0, mi) and the mode 2 docs 
(n, £, to) = (ri2, 0, 7712) . We also define the tilde matrix 



A Q B 
-B -Aq 



(57) 



Next, we consider the j ' ^ part of the free Hamilto- 
nian (|25p. In the matrix representation, H J under the 
two-mode approximation is written as 



\z),H f Zy 



(58) 
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where 



Z v = 



z\,= 



Uy 



Ay 



B*, = 



0,2' 



\4>J 

a[, a 2 i ay ay 

Ay By 
By Ay t 



2X(1',2 / ) 
2K(2',V) 





£ 2 /-M + 4i<r(2',2 / ) 



(59) 

(60) 
(61) 

(62) 
(63) 



Here, the quantum numbers (n,£,m) of the mode 1' and 
the mode 2' are {n{,j',m'{) and (n' 2 , —j' ,m' 2 ), respec- 
tively. The tilde matrix is defined by 



Uy 



Ay By 

-By -Ay 



(64) 



B. Small coupling expansion 

In this subsection, we define the dimensionlcss param- 
eter A as 



A = d ho (gN c ), 



and make an assumption of 



|A| « 1, 



(65) 



(66) 



where o?ho = (a^a 2 ) 1 / 3 M is the constant determined by 
the typical scales of the system, and expand the func- 
tion f(r, z) (see Eq. iJTIJ) and chemical potential /x with 
respect to this small parameter A. We refer to this expan- 
sion as the small coupling expansion. Here we perform 
this expansion up to the first order in A. Note that this 
expansion is applicable not only for repulsive interaction 
(A > 0), but also for attractive one (A < 0). Let us intro- 
duce the notations of quantities of the zeroth order in A, 
namely, a real function fo(r, z) and a chemical potential 
fj®, which satisfy 



[K K + V(r,z)}f (r,z)=nofa(r,z). 



(67) 



The function fo(r,z) and the quantity /io correspond to 
the eigenfunction wo,o,o( r 7 z) and the eigenvalue £o,o,o, re- 
spectively, i.e., 



fo(r, z) = V27ruo,o,o(r-,2;) 

= V2^C 0fifi e-^ r2 (a ± r) K e-^ z2 , (68) 

(69) 



Note that the function fo(r, z) is normalized as 

rdrdzf$(r,z) = 1. (70) 



Using the function fa(r, z) we can expand the real func- 
tion /(r, z) up to the first order in A: 

f(r,z)~G{f (r,z) + \h(r,z)} , (71) 

where fi(r, z) is a real function and normalized, 

rdrdzf 1 (r,z) = l, (72) 



and the parameter G is determined by the normalization 
condition (fLl)) as 

G = (1 + 2X Pl )-i , (73) 
where the coefficient pi is defined as 

Pi = Jrdrdzf (r,z)fi{r,z), (74) 

which leads to 

f(r, z) ~ / (r, z) + A{/i(r, z) - Pl f (r, z)} . (75) 

We also need to expand the chemical potential /x up to 
the first order in A: 



(76) 



Mo = £o : o,o = ^i( K + 1) + ~2 ■ 



We determine \i\ by substituting Eq. (|75|l and Eq. Ij76(l 
into the GP equation l |T3j l. 



{K K +V(r, z)-Mo}/i(r, z)-^f Q {r, z)+_— / Q 3 (r, z) = . 

(77) 

Here we expand the function /i (r, z) with the cigenfunc- 
tions u n< o,m(r, z), 



fx(r,z)= ^ Sn,mU n ,0,m(r, z) , (78) 
n,m— 

where the expansion coefficients s„ !m are written as 

s n ,m = 27T / ' rdr dz u n ,o, m {r, z) /i(r, z) . (79) 

It is notable to see that 

s ,o = V^TTPl • (80) 
Substitute Eq. JTSJ to Eq. lff7)l. one obtains 

OO 

W n) o,m(^) z ) - Mi/o(^ 2;) 

n,m— 

+ ^^-/o(r,«) = 0. (81) 



271-d 



ho 
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Multiplying the function fo(r, z) to Eq. I|81[) and inte- 
grating over the space coordinates, we obtain 

*-5^/'**>«'-'>-£ c !$ (82 » 

where C is defined as 
C = 



1 aj_a 2 



2?r 5 2 2k + 2 
Thus fi(r, z) and /^i are determined. 



(83) 



C. The condition for the existence of the complex 
modes 

In this subsection, we calculate the eigenvalues of the 
matrices Tio in Eq. (|57Jl and Hj> in Eq. {f&^i under the 
small coupling expansion. To simplify the expression, 
we parameterize the elements of the matrices Aj and Bj 
C? = 0,j')as 



4»,n 


^3 ~^ ^j: a 1 


(84) 


A S,22 


= aj — rj a , 


(85) 




= Aj t 2i = Aj, , 


(86) 


B iM 


= b 3 + r 3,b , 


(87) 


B 3,22 


= b 3 - r 3,b , 


(88) 


B j,12 


= Bjfli = Aj, 6 . 


(89) 



The eigenvalues Xj °f t ne matrix Tlj are characterized by 
the equation 

det[Hj- X j] =0. (90) 

Then we obtain 

Xi = K 2 + + A D - ( b2 3 + r h + A h) ± 

(91) 

7 = (a>jrj, a - h r 3\bf + (aj Aj, a - bjAj, b ) 2 

- (rj, a Aj, b - rj, 6 Aj, a ) 2 . (92) 
We discuss the two cases, j = and j ^ 0, separately. 

i. j = case 

In the case of j = 0, the elements of the matrix Aq are 
written as 

a = -{2ci)j_(ni + ra 2 ) + w z (mi + m 2 )} + 0(A) (93) 



= a + O(A), 



(94) 



ro, a = -{2w\_{ni - n 2 ) + uj z (m 1 - m 2 )} + O(A) (95) 



= f , a + O(A) . 

It is easy to see that if |A| <C 1, then 

X 2 , = (ao ± r , a ) 2 + O(A) > . 



(96) 
(97) 



^From this relation, we conclude that no complex mode 
appears for j — and sufficiently small A. 



2. j 7^ case 
In the case of j ^ 0, it is straightforward to see 

X 2 ? = (r jl , a ± y j4-A^ b ) 2 (98) 

for j' > 1. So the sign of a 2 / — A 2 , b determines whether 
the eigenvalues are real or not. Note that it has the term 
of ey +£2', and when we search complex eigenvalues, it is 
natural to choose the mode 1' as (0, j',0) and the mode 
2' as (0, — j',0), because then a 2 , reaches its minimum. 
In the case of j' > k, sy and e 2 > are written as 



ey = U ± (j' + K+l) + 

e 2 > = u>±(j' -« + !) + 



(99) 
(100) 



Since the matrix element a* is written as 



X/j, 1 +2K(l',l')+2K(2',2'), (101) 



the term a 2 , — A 2 , & is evaluated as 



, 2 ,-A 2 ^ b = c 2 (/- K ) 2 + O(A)>0. 



(102) 



So we can see that in the case of j' > k no complex mode 
exists. 

In the case of < j' < n, ey and e 2 i are 



ey =CJ X (f + K + l) + y , 

e v =u±(-j' + K + l) + Y' 
Then the matrix element aj' becomes 

aj , = -A/ii + 2K(1', 1') + 2K(2', 2') 
where K (!', 1') and K(2', 2') are written as 



K(l\l>) = —\ 



2tt C (2k + j 



Ml 



rfho 2 2J'(k +/)!«! ' 
, , „ 2tt C (2k-/)! 

2', 2' - — A— -— jr- ^— 

dho 2 2 ■? (k — 
The parameter Aj' jb is given as 

A,, b = ^AO- ^ 



d ho kV(k + /)K«-J')J 



(103) 
(104) 

(105) 

(106) 
(107) 

(108) 



It is important to see that aj' and Aj',6 have no term of 
the zeroth order in the parameter A (cf. Eqs. I|93|l and 
(11 IB ) . Thus we must evaluate all terms in a 2 / — A 2 , b up 
to the second order in A. Then, a 2 , — A 2 , b is evaluated 



a,-A 2 =A 2 M 



2tt 
2tt 



C(S + T) 



(109) 
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3 


4 


5 




1 














2 




X 










3 




X 


X 








A 
4 




X 


X 








5 




X 


X 


X 
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X 


X 


X 





















TABLE I: The existence of the complex modes. The condition 
11121 1 are examined for various values of (k, j). The marks are 
the positions where the complex modes arise. This result is 
valid for the repulsive interaction (A > 0), as well as the 
attractive one (A < 0), under the situation of |A| 1. 



where S and T arc defined as 

{2K + j')\ , (2K-j')\ 



S 
T = 



(110) 
(111) 



Therefore, the condition for the existence of the complex 



modes, af/ — Af/ b < 0, now amounts to 



5' 



(«!) 



h2 



T. 



(112) 



We apply the condition i|112|) to each winding number 
and find the appearance of the complex modes (Table 
QJ. In the case of the winding number k = 1, no com- 
plex mode exists. On the other hand, in the case of the 
winding number k > 2 the complex modes always exist. 
Particularly, we find that in the case of the winding num- 
ber k = 2 the complex modes appear as a pair, while in 
case of k — 3 the complex modes arise as two pairs. So 
our results in the case of the winding number k — 2 and 3 
are consistent with those by Pu et al. 01 an d Mottonen 
et al. ^lj. I n the case k = 4, our result is consistent with 
tha t by Kawaguchi and Ohmi Remark that Pu et 

al. [l(j performs their analysis also for the attractive in- 
teraction (A < 0), while the other authors do only for the 
repulsive interaction (A > 0) 0, R^l ■ O ur result is also 
consistent with that by Pu et al. [l(| for A < 0. 

The condition (jTT2j) can be applicable only in the small 
coupling constant region. So the existence of new com- 
plex modes which may arise in the large coupling con- 
stant region, such as one can find them in Fig. 2 in 
Rcf. 12], is outside the scope of this paper. 



D. Proof of the existence of the complex modes for 

k > 2 

One of the most important consequence of the analytic 
result 1112fl is the following statement: The j' = 2 modes 



are always complex ones for k > 2. The proof is given 
below. 

For j' = 2, the inequality l|112f) is reduced to the fol- 
lowing coupled inequalities 



1 < 



(2k + 2)(2k+1) 4k(k-1) 



k(k — 1) 



4(k + 2)(k + 1) (2k)(2k-1) V (k + 2)(k + 1) ' 

(113) 



(2k + 2)(2k+1) Ak(k-I) 



k(k — 1) 



4(k + 2)(k+1) (2k)(2k-1) V(k + 2)(k + 1) 



< 1 



(114) 



First, let us prove the inequality (|113fl . Remark that 
the inequality 



(2k + 2)(2k+1) 4k(k-1) 



4(k + 2)(« + 1) (2k)(2k-1) 
is reduced to 

k(k + 1)(4k 2 -2k-5) >0. 



- 1 > 



(115) 



(116) 



The largest root for the left-hand side is k = 
K\ ~ 1.396. Combining this result with the fact 
y/ {k(k — 1)}/{(k + 2)(k + 1)} > for k > Ki, we can 
find the inequality (|113|l is satisfied for k > k\. 

Next, let us prove the inequality 1114(1 . The inequality 
(|114fl is equivalent to 



k 2 (k + 1)(8k 3 + 52k 2 - 53k - 25) > . 



(117) 



The largest root for the left-hand side is k = K2 — 1.199. 
^From the same discussion on the inequality (|113|) . we 
find that the inequality (|114|1 is satisfied by n > K2- 

^From our discussion above, the roots k,\ and satisfy 
the inequality 1 < Ki, k% < 2. This means that if n > 2, k 
satisfies the inequalities (111 3|> and (|114|) simultaneously, 
then the statement is proven. 



V. SUMMARY 

We derived the analytic expression of the condition 
for the existence of the complex modes when the con- 
densate has a highly quantized vortex, using the method 
developed by Rossignoli and Kowalski 0] for the small 
coupling constant, under the two- mode approximation. 
Our results agree with those by Pu et al. (k — 2, 3) 
[Toj . Mottonen et al. (k = 2) [ll| and Kawaguchi and 
Ohmi (k = 4) 12]. It is emphasized that the formula 
derived here is applicable for an arbitrary winding num- 
ber, and one can check whether complex modes exist or 
not for an arbitrary highly quantized vortex in the con- 
densate. We have also proven that in the case of k > 2, 
the complex modes always exist. This means that the 
arbitrarily highly quantized vortex may be dynamically 
unstable. Moreover, it is shown that imaginary parts of 
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the complex eigenvalues rise linearly against the coupling 
constant. 

The validity of the two-mode approximation is cer- 
tainly not clear. But the numerical calculation for the 
case of k = 2 and A > revealed that the behavior of 
the imaginary part in the small coupling constant region 
is not significantly affected by the number of the modes 
included to the calculation |10|. This fact suggests that 
the two-mode approximation does not modify behaviors 
of complex modes in the small coupling constant region 
essentially and that the analysis in this paper is reliable. 

Finally, we comment that the physical interpretation 
of the complex modes is very difficult and is not settled 
yet. Pu et al. employs the classical interpretation on the 
complex modes, with the idea that the complex modes 
cause the dynamical instability, that is, the decay of the 
initial configuration of the condensate [ljj. The argu- 
ment does not touch the aspect of quantum theory. We 
recently proposed a quantum field theoretical treatment 



associated with the complex modes and sought a new 
interpretation |19| . 
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